Dissipative quantum chaos: transition from wave packet collapse to explosion 
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Using the quantum trajectòries approach we study the quantum dynamics of a dissipative chaotic 
system described by the Zaslavsky map. For strong dissipation the quantum wave function in the 
phase space collapses onto a compact packet which follows classical chaotic dynamics and whose area 
is proportional to the Planck constant. At weak dissipation the exponential instability of quantum 
dynamics on the Ehrenfest time scale dominates and leads to wave packet explosion. The transition 
from collapse to explosion takes place when the dissipation time scale exceeds the Ehrenfest time. 
For integrable nonlinear dynamics the explosion practically disappears leaving place to collapse. 
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Technological progress leads to the investigation of 
physical phenomena at smaller and smaller scales, where 
both quantum and dissipative effects play a very impor- 
tant role. At present, general theoretical concepts for 
the description of quantum dissipative systems are well 
developed and established [1, 2]. A major tool is the 
master equation, that governs the evolution of the den- 
sity matrix [3]. For the simplest dynamics this equation 
can be solved exactly. However, for complex nonlinear 
systems analytical solution is absent and even numerical 
simulations become very difhcult. Indeed, for a system 
whose Hilbert space has dimension N, one has to store 
and evolve a density matrix of size N x N. In spite 
of these limitations, numerical simulations of the mas- 
ter equation allowed to perform the first studies of the 
quantum dynamics of classically chaotic dissipative sys- 
tems showing a quantum strange attractor [4]. 

Quantum trajectòries are a very convenient tool to sim- 
ulate dissipative systems [5, 6]. Instead of direct solu- 
tion of the master equation, quantum trajectòries allow 
us to store only a stochastically evolving state vector of 
size N. By averaging over many runs we get the same 
probabilities (within statistical errors) as the ones ob- 
tained through the density matrix directly. Besides their 
practical convenience, quantum trajectòries also provide 
a good illustration of individual experimental runs [7]. 
Indeed, modern experiments often enable us to address 
a single quantum system evolving under the unavoidable 
influence of the environment. 

It is known that for linear systems dissipation leads to 
wave packet localization [8] . Numerical results as well as 
theoretical arguments indicated that localization can oc- 
cur also in nonlinear systems [9, 10]. On the other side, in 
absence of dissipation it is known that the instability of 
classical dynamics leads to exponentially fast spreading 
of the quantum wave packet on the logarithmically short 



Ehrenfest time scale íg ~ |ln?i|/À [11, 12]. Here A de- 
notes the Lyapunov exponent which gives the rate of ex- 
ponential instability of classical chaotic motion, and h is 
the dimcnsionless effective Planck constant of the system. 
In this paper we show that for the dissipative quantum 
chaos there exist two regimes: one with the wave packet 
explosion (delocalization) induced by chaotic dynamics 
and another with the wave packet collapse (localization) 
caused by dissipation. We argue that the transition (or 
crossover) from collapse to explosion takes place when the 
dissipation time I/7 becomes larger than the Ehrenfest 
time íe (7 is the dissipation rate). 

We investigate the quantum evolution of a kicked sys- 
tem subjected to a dissipative friction force. Assuming 
the Markov approximation, we can write a master equa- 
tion in the Lindblad form [3] : 



P : 



(1) 



where p is the density operator, { , } denotes the anticom- 
mutator, are the Lindblad operators, which model the 
effects of the environment, and H is the Hamiltonian of 
the system. We consider a kicked system, described by 
the Hamiltonian 
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where T is the kicking period and the operators x and 
íi = —i(d/dx) come from quantizing the classical vari- 
ables x G [0, 27r[ and n £ (— 00, +00). This Hamiltonian 
corresponds to the kicked rotator [13], a paradigmàtic 
model in the fields of nonlinear dynamics and quantum 
chaos. This model is also on the focus of experimen- 
tal investigations with cold àtoms in optical lattices [14]. 
We assume that dissipation is described by the lowering 
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operators 
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with n = 0, 1, ... eigenvalues of the operator n. At the 
classical limit, the evolution of the system in one period 
is described by the Zaslavsky map [15] 



n t+ i = (1 - 7)n t + ksinxt, 
x t +i = x t + Tnt+i, 



(4) 



where the discrete time t is measured in number of kicks 
and 1 — 7 = exp(— g 2 ). This map describes a friction force 
proportional to velocity. We have < 7 < 1; the limiting 
cases 7 = and 7 = 1 correspond to Hamiltonian evolu- 
tion and overdamped case, respectively. Introducing the 
rescaled momentum variable p — Tn, we can see that 
the classical dynamics depends only on the parameters 
K = kT and 7. Since [x,p] — [x, Tn] — iT, the effective 
Planck constant ish = T. The classical limit corresponds 
to H — > 0, while keeping K — hk — const. 

The first two terms of Eq. (1) can be regarded as the 
evolution governed by an effective non-Hcrmitian Hamil- 
tonian, H cS = H s + iW, with W = -l^E^iV In 
turn, the last term is responsible for the so-called quan- 
tum jumps. Taking an initial state |</>(to)), the jump 
probabilities dp^ in an infinitesimal time dr are defined 
by dpp_ = (0(ro)|ït£/i|0(To))dT, and the new states af- 
tcr the jumps by |^) = L^{t )) /\\L^{t ))\\. With 
probability dp^ a jump occurs and the system is left in 
the state . With probability 1 — there are 

no jumps and the system evolves according to the effec- 
tive Hamiltonian ií c ff. In this case we end up with the 
state = (1 - iH eS dt/h)\cf)(to))/y/Í - £ fe d Pk . We 
note that the normalization is included also in this case 
because the evolution governed by H r is non-Hermitian. 
To simulate numerically the above described jump pic- 
ture we use the so-called Monte Cario wave function ap- 
proach [7] . The changing state of a single open quantum 
system is represented directly by a stochastically evolving 
quantum wave function, as for a single run of a laboratory 
experiment. We say that a single evolution is a quantum 
trajectory. 

We focus first on the chaotic regime for the kicked rota- 
tor dynamics. Therefore, we consider K = 7, correspond- 
ing to a positive Lyapunov exponent A « ln(if/2) = 1.25. 
The localization-delocalization transition is clearly illus- 
trated in the two top panels of Fig. 1. They show the 
Husimi function [16] corresponding to a single quantum 
trajectory evolution, computed after t — 300 kicks. In 
both cases the initial wave packet is a Gaussian state 
with equal uncertainties Aa; = Ap = y/H/2. We can see 
that for strong dissipation (7 = 0.5) the wave function 
of a single quantum trajectory at t — 300 is localized in 
the phase space (top left panel in Fig. 1). In contrast, 
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FIG. 1: (color online) Top: Husimi functions in phase space 
for a single quantum trajectory taken after t = 300 kicks, 
at K = 7, % = 0.012, 7 = 0.5 (left) and 7 = 0.01 (right). 
Here the coordinates x (horizontal axis) and p (vertical axis) 
vary in the intervals: < x < 2n, —25 < p < 25 (left) and 
— 100 < p < 50 (right); the width of the p-interval is the same 
in both cases for comparison purposes. The initial Gaussian 
wave packet is located at ({x),{p}) = (57r/4, 0). The color is 
proportional to density: blue for zero and red for maximum. 
Bottom: quantum Poincaré section (left), obtained from av- 
erage quantum x,p vàlues for the case of top left panel and 
its classical counterpart (right). In these panels < x < 2n 
and -15 < p < 15. 



the case of weak dissipation (7 = 0.01, top right) is char- 
acterized by wave packet delocalization. Since for strong 
dissipation the wave packet is localized in phase space, 
it makes sense to draw a quantum Poincaré section by 
printing the expectation vàlues (x) and (p) at each map 
step. The quantum Poincaré section is shown in Fig. 1 
(bottom left) and is characterized by the appearance of 
a strange attractor. A very similar strange attractor is 
obtained also from the classical Poincaré section corre- 
sponding to the Zaslavsky map (4) (see Fig. 1 bottom 
right). We also note that the Husimi function obtained 
in the case of weak dissipation exhibits a spreading of the 
quantum wave packet over the strange attractor. Also in 
the strongly dissipative regime the localized wave packet 
is stretched along the direction of the attractor. 

A further confirmation of the good agreement between 
the classical and quantum dissipative evolutions is ob- 
tained by computing the function / = {p)t+i — (1 ~l){p)t 
[17]. From the classical map (4) we expect f(x) = 
K since. The comparison between f((x)) and the func- 
tion f(x), shown in Fig. 2, indicates that the Zaslavsky 
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FIG. 2: Comparison between f({x)) computed from the quan- 
tum dynamics (circles) and the kick function f(x) = ifsina; 
(sòlid curve). Parameter vàlues and initial conditions are as 
in Fig. 1 (left top and bottom panels). Error bars represent 
the quantum uncertainty in f({x}). 



map provides a good description of the quantum wave 
packet motion. Indeed the points /((x)) are concentrated 
around the curve f(x), with dispersion proportional to 
y/h (see Fig. 4 below). Therefore a quantum trajectory 
exhibits the same important features of a classical tra- 
jectory, including the exponential instability, with rate 
given by the classical Lyapunov exponent. A significant 
difference between quantum and classical trajectòries is 
the presence of quantum noise [18]. Therefore, a more 
precise identification can be done between quantum evo- 
lution and noisy classical evolution, the noise amplitude 
being Ax ~ Ap oc y/h. It is interesting to note that the 
chaotic behavior of classical systems can be reproduced 
also by non-dissipative continuously measured quantum 
systems [19, 20, 21]. 

The wave packet dispersion is measured by a t = 
yj (Ax)'j + (Ap)'j. This quantity is evaluated, for weak 
and strong dissipation, in Fig. 3 (left panel), using the 
same parameter vàlues and initial conditions as in Fig. 1. 
In both cases there are strong fluctuations, which we 
smooth down by computing the cumulative average a t = 
j ^2j—i Cj- The convergence of the time averaged quan- 
tity ef t to a limiting value is clear. It is also evident 
that the wave packet spreading is much stronger at weak 
than at strong dissipation. We would like to stress that 
the same limiting value of the average dispersion a is 
obtained for any quantum trajectory, independently of 
the initial condition. This is demonstrated in the right 
panel of Fig. 3, where we compare at for two completely 
different initial conditions: a Gaussian wave packet and 
an eigenstate of the operator x, that is, \x) = \xq). In 
the latter case there is a complete delocalization along p 
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FIG. 3: Left: Dispersion a (fluctuating curves) and cumula- 
tive average a (smooth curves) as a function of the number 
of kicks t, for the cases considered in Fig. 1, with 7 = 0.01 
(dashed curves) and 7 = 0.5 (full curves). Right: a versus 
time t at 7 = 0.5, starting from an initial Gaussian wave 
packet (full curve, also shown in the left panel) and a coordi- 
nate eigenstate \x) = \%o), with xo = ir (dashed curve). 



(limited only by the size of the basis considered in our nu- 
merical simulations) and dissipation leads to the collapsc 
of the wave packet, which eventually becomes localized 
in phase space. 

The average dispersion a of the wave packet as a func- 
tion on the dissipation strength 7 is shown in Fig. 4, 
for a few vàlues of the effective Planck constant h, with 
0.012 < h < 0.33. The localization-delocalization tran- 
sition can be seen for all h vàlues. In Fig. 4 inset we 
consider the sealed dispersion a s = W/ y/h. At strong dis- 
sipation all curves collapse, while at weak dissipation the 
sealing a oc y/h is not fulfilled. Our numerical results can 
be explained as follows. Due to the exponential instabil- 
ity of chaotic dynamics the wave packet spreads exponen- 
tially and for times shorter than the Ehrcnfcst time we 
have a t ~ y/hexp(Xt). The dissipation localizes the wave 
packet on a time scale of the order of I/7. Therefore, for 
I/7 <C íe, we obtain a ~ y/Kexp(X/j) <C 1. In contrast, 
for 1 /7 > íe the chaotic wave packet explosion dominates 
over dissipation and we have complete delocalization over 
the angle variable. In addition, in this case, the wave 
packet spreads algebraically due to diffusion for t > íg. 
For t > t E we have a t ~ y/D{K)t, D(K) w K 2 /2 being 
the diffusion coefficient. This regime continues up to the 
dissipation time I/7, so that a ~ y/D(K)/j. Accord- 
ing to the above argument, the transition from collapsc 
to explosion, which we wish to call Ehrenfest explosion, 
takes place at 



t E ~ |lnft|/A~ I/7 



(5) 



Our numerical data at moderate vàlues of h > 0.01 indi- 
cate a smooth transition (crossover). However, we can- 
not exclude from our data that in the limit h — > the 
transition becomes sharp. Since the dependence on h is 
only logarithmic it is difhcult to check numcrically the 
above relation. However, it is compatible with our data 
obtained for h > 0.01. First of all, in the localized regime 



4 




0.2 0.4 0.6 0.8 1 



Y 

FIG. 4: Average dispersion a as a function of 7, for K = 7, 
fi = 0.33 (circles), 0.11 (squares), 0.036 (diamonds) and 0.012 
(triangles). Stars show the same quantity for the integrable 
case K = 0.7, at % = 0.012. Inset: scaled dispersion a s = 
a/Vh versus 7, with same meaning of symbols as in the main 
figure. 

7Í£ > 1 the scaling law a oc \fh is satisfied. Moreover, it 
is satisfied down to smaller and smaller 7 vàlues when h 
is reduced. Therefore, even for infinitesimal dissipation 
strengths the quantum wave packet is eventually local- 
ized when h — > 0: we have lim^o a = 0. In contrast, in 
the Hamiltonian case (7 = 0) linift_+o<5 : = 00. This re- 
sult underlines the importance of a (dissipative) environ- 
ment in driving the quantum to classical transition: only 
for open quantum systems the classical concept of tra- 
jectory is meaningful for arbitrarily long times. On the 
contrary, for Hamiltonian systems a description based on 
wave packet trajectòries is possible only up to the Ehren- 
fest time scale. 

We would like to emphasize the role played by chaotic 
motion. For this purpose, in Fig. 4 we also show a as 
a function of 7 in the integrable regime at K = 0.7, for 
h = 0.012. In this case the wave packet dispersion is 
much smaller than in the chaotic regime: the Ehrenfest 
time scale is algebraic and not logarithmic in h. Thus, a 
much weaker dissipation is sufficient to localize the wave 
function in the case of integrable dynamics. This can be 
clearly seen from our numerical data shown in Fig. 4. 

Finally, we would like to point out that the transi- 
tion described here could be observed by means of Bose- 
Einstein condensates in optical lattices. A first experi- 
mental implementation of the kicked rotator model us- 
ing a Bose-Einstein condensate has been recently re- 
ported [22]. Dissipative cooling techniques are possible 
in these systems. Moreover, images of atòmic clouds can 
be taken, thus measuring their dispersion. Also the mea- 
sured condensate positions should give a clean kick func- 



tion (like in Fig. 2) in the case of collapse and random 
scattered points in the case of cxplosion. Such experi- 
ments would give important information not only on the 
interplay between chaos and dissipation but also on the 
stability of the condensate [23] under the joint effects of 
chaotic dynamics and dissipation. 
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